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We report the results of a numerical simulation of a lattice formulation of the two-dimensional 
= (2,2) super Yang-Mills theory proposed by Suzuki and Taniguchi [jl|]. We measure the 1- 
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a specific value gives a supersymmetric continuum limit. Our results are consistent with this 
scenario although conclusive results on the restoration of supersymmetry have not been obtained. 
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1. Introduction 

A lattice formulation of the supersymmetric gauge theory is important to understand the non- 
perturbative aspects. Recently some formulations are proposed by several authors Most of 
the models use the fact that a nilpotent part of the supersymmetry can be kept on a lattice for 
N >2 cases using the topological twist and the relations among some of them are now becoming 
transparent ^, ||]. There is also an attempt to keep all the supersymmetry |7|]. Another 
approach is a model without exact supersymmetry [|l|] which uses the fact that in two-dimensional 
case, because of the super-renormalizability, there are only a few fine-tuning parameters to obtain 
the supersymmetric continuum limit. One of the merits of lattice formulation is to enable one to 
perform numerical simulations. In two-dimensional case, some results are known for the super 
Yang-Mills theory with the topological twist ||8|, 

In this talk, we report the result of a numerical simulation of the model without exact super- 
symmetry proposed in [[I]]. We measure some 1 -point and 2-point functions. 

2. Model and Algorithm 

The target theory in the continuum is the 2-dimensional = (2, 2) supersymmetric Yang-Mills 
theory which is obtained by a dimensional reduction from the 4-dimensional A'^ = 1 super Yang- 
Mills. The lattice action is defined as a lattice version of a 4-dimensional action onLxLxlxl 
lattice together with a scalar mass counter term 5counter:^ 

S = Sg + Sf + '^counter' 

(2.1) 

The bosonic part is a plaquette action 

5G[t/] = :r^£ £Retr{l-P(x,M,A^)}, (2.2) 

P{x,M,N) = U{x,M)Uix + aM,N)Uix + aN,M)-^U{x,N)-\ (2.3) 

where we use SU{Nc) gauge link variables, U{x,}Jl) = &xp{agA"^{x)T'^), (/i = 0, 1) and compact 
scalar fields, U{x,2) = exp(a^(p"(x)r") and U{x,'i) = exp{ag(j)" {x)T"). The coupling constant g 
is related to j8 through j3 = INc/cP'g^- The fermion action consists of the Wilson-Dirac operator 

1 3 

5f[[/,A] = -a^ £tr{ A (x)CD^ A (x)}, Z)^ = - £ {rM(Vl, + Vm) -«V1,Vm}, (2.4) 
xer ^ M=o 

with covariant differences for the adjoint representation 

VmA(;c) = ^{U{x,M)X{x + aM)U{x,My^ -X{x)] (2.5) 
and its adjoint . The counter term is 

5counter[t/] = "^A^. £ (tr{[/ (x, 3) + [/ (x, 3) - 2} + tr{[/(x, 2) + C/(x, 2)" ^ -2}) , (2.6) 

-ver 



'For the notational details, see 
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where = 0.65948255(8). The counter term is intended to cancel the radiative corrections to the 
scalar mass term. Other corrections which might appear in the effective action are suppressed in 
the continuum limit because of the super-renormalizability of this theory. All possible divergences 
in the sub-diagrams in the perturbative expansion are suppressed. It should be noted, however, that 
it does not guarantee the supersymmetry of composite operators. 

In our numerical simulation, we use quenched gauge configurations generated by using the 
Hybrid Monte Carlo algorithm. The fermion contribution is introduced as a reweighting by the 
Pfaffian. In the continuum limit, the model has a real and positive Pfaffian. In fact the direct calcu- 
lation for some sample configurations shows that the Pfaffian is real and positive in our parameter 
region. Therefore we use a positive square root of the determinant which numerical cost is much 
less expensive. We set a bare fermion mass m = 0. The lattice size is 8 x 8 for 1 -point functions 
and 12 x 12 for 2-point functions. Since the coupling is j8 = 2Nc/a^g^, the continuum limit is the 
j5 oo limit. We set 3 < j3 < 40. The gauge group is SU (2). 

We summarize the parameters and the numbers of independent configurations in table |l|. 
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Table 1: The numbers of configurations for each parameter set on 8 x 8 lattice. 



3. One-point functions 

It is of our interest to measure the vacuum expectation value of a supercharge-exact operator, 
{Qff), because it must vanish in the supersymmetric continuum limit. It should be noted, however, 
that can be non-zero (a finite renormalization), depending on the definition of the composite 
operator Qd' which does not necessarily preserve the supersymmetry. 

We make use of a scalar part of the topological twisted supercharges Q and observe Q-exact 1- 
point functions used in [||, Since we have no exact supersymmetry at finite lattice spacings, first 
we write down the continuum relations and then discretize them. We define the scalar supercharge 
Q in the continuum as follows: 

QA"^ = Vt,, Qwt,=iD^^'^, ef" = 0, (3.1) 

Q^'" = r]^ e Y = -^^/..cf , QX" = iFSx . (3.2) 

Here we introduce scalar fields = (p + i(^ and (j)' = (p — icj) . Fermions in the twisted basis y/^, rj 
and X are given by liner combinations of the components of A: {ri/2,x,Yo,¥iV = T^- 
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We use the following three ^,'s: 

^1 = -^§r]%/.c0 V, ^2 = -2a"foi , ^3 = -^r^D^'^'"- (3.3) 

We divide Q^,- into two parts Q^i = Fi + Bi, where Fj contains fermions and Bi is made only from 
bosons: 

= -^gfatcr^'^" , F2 = ix"{DoW?-DyWS) , ^3 = ^r^D^^ + ^gfabc'^'^Xj/'^rn , 

Br = -^g' {fabc^"^''y, B2 = 1{FS,)\ B, = ^D.rD,-^'". (3.4) 

In the continuum, the kinetic terms and Yukawa interactions of the twisted fermions are contained 
in the Dirac operator. Therefore, to find an appropriate discretization of these terms, we simply 
replace the continuum Dirac operator with the Wilson-Dirac operator. For example, we replace 

iVoDo^ Vo {{T-'f{CDw)T-')^^^^ |. (3.5) 

All dimensionful observables are measured in a unit of the dimensionful coupling g. 

First, we present the result of the case. Figure |l] shows that each of the bosonic and 
fermionic parts is divergent in the continuum limit. The sum {Qff\) stays finite after the reweight- 
ing by the Pfaffian, while the quenched result diverges (Fig. In this theory, the cancellation of 
divergences in {QG\) is achieved by a balance between bosons' and fermions' degrees of freedom. 
Our reweighted result is consistent with this fact and the effect of dynamical fermions appears to 
be properly included by the reweighting. As already noted, even in the supersymmetric contin- 
uum limit, (2^1 ) can be non-zero due to a finite renormalization. The dependence of {QG\) 
is summarized in Figure ||. To estimate the effect of the finite renormalization, we determined the 
values of which provide {Q&x) = 0. They are <^ = 1.047(51) at jS = 8 and = 1.006(77) 
at j8 = 20. Almost no j8 -dependence is observed. These values are significantly different from 
'ta = 0.65948255 calculated in the continuum limit [|l|] and suggest that the effect of the finite 
renormalization is certainly not negligible. 

The cancellation of divergences is also realized both in Qff2 and Q^j, cases (Fig. The 
"^-dependence, however, is not manifest. Note that the plots have much more errors than that of 
{QG\) and this would imply that the "^-dependence is smeared. The difference between behavior 
of QG\ and that of Q62 and Q6-i could be accounted as a result of the difference of the divergence 
of each operators. B\ and F\ have logarithmic divergences, while B2, B3, F2 and F3 have quadratic 
divergences. 



4. Two-point functions 

Next we measure quantities with which the restoration of supersymmetry is expected to be 
observed transparently. A supersymmetric Ward-Takahashi identity indicates that the following 
2-point functions should have the same functional form:^ 

B = 2i{h,ix)My)) , F = {tiir.Ysir + ir5r)ri^)j7'\y)}) , (4.1) 

^The former function B exhibits a power law behavior in the continuum theory llQ]. 
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Figure 1: The a^O limit of Bi and Fu the left is 
is ^ = 0.65948255. The lattice size is 8 x 8. 
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Figure 2: The a ^ limit of {Q0\), the left is reweighted and the right is quenched. The lattice size is 



where the bosonic currents are 

i^{x)=r'Yp.v{^) , j5^{x)=r'Ynr5v"i^)+2i{rdn(p"ix)-(p"d^rix)}, (4.2) 

and the fermionic current is 

W = TYu { ^F^aOpo - iypDp ir + - i8fabc<p'^'y5 } W . (4.3) 

Figure |5| shows a typical result of the 2-point functions. The number of the configurations we used 
is 101 and the bare fermion mass is 0. According to the scenario, we expect that a suitable choice 
of should give the supersymmetric result, i.e., the identical spectra, while the other choices of 
should not. Unfortunately, errors in the plot are too large to analyze the spectra although this is a 
result with the quenched approximation. The point here is that we cannot distinguish the difference 
of the counter term. 

5. Conclusion 

We observed 1 -point functions and 2-point functions in a lattice formulation of the two- 
dimensional N = (2,2) super Yang-Mills theory. In our scenario, only the counter term coefficient 
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Figure 4: The dependence of {QO2) (left) and {QO3) (right). These are reweighted and the lattice size is 



should be finely tuned. The 2-point functions have large errors that we cannot compare the spec- 
tra associated with a bosonic current and a fermionic current. The 1 -point functions we observed 
are finite in the continuum limit because of the fermion loop effect. The result of a less divergent 
1 -point function depends on and is consistent with our scenario. To obtain the conclusive result 
from this dependence, i.e., whether the scenario actually works or not, we need the renormaliza- 
tion factor for the 1 -point function. Although the current result is not quite promising, we have 
some possible ways to improve. A UV-filtered reweighting will help to reduce the errors after the 
reweighting. The HMC algorithm with dynamical fermions is another option. The result of 2-point 
functions suggests that the fermion or the scalar (or both) is rather far from massless so that a neg- 
ative bare mass of the fermion which reduces the physical mass may improve the sensitivity on the 
counter term. 
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